Gutzwiller Projected wavefunctions in the fermonic theory of S=l spin chains 



Zheng-Xin Liu, 1 Yi Zhou, 2 Hong-Hao Tu, 3 Xiao-Gang Wen, 4 ' 1 and Tai-Kai Ng 5 

1 Institute for Advanced Study, Tsinghua University, Beijing, 100084, P- R- China 
2 Department of Physics, Zhejiang University, Hangzhou, 310027, P.R. China 
3 Max- Planck- Institut fiir Quantenoptik, Hans-Kopfermann-Str. 1, 85748 Garching, Germany 
4 Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA 
5 Department of Physics, Hong Kong University of Science and Technology, Clear Water Bay Road, Kowloon, Hong Kon> 

We study in this paper a series of Gutzwiller Projected wavefunctions for S — 1 spin chains 
obtained from a fermionic mean-field theory for general S > 1/2 spin systems [Phys. Rev. B 81, 
224417] applied to the bilinear-biquadratic (J-K) model. The free-fermion mean field states before 
the projection are ID paring states. By comparing the energies and correlation functions of the pro- 
jected pairing states with those obtained from known results, we show that the optimized Gutzwiller 
projected wavefunctions are very good trial ground state wavefunctions for the antiferromagnetic 
bilinear-biquadratic model in the regime K < J, (J > 0). We find that different topological phases 
of the free-fermion paring states correspond to different spin phases: the weak pairing (topologically 
non-trivial) state gives rise to the Haldane phase, whereas the strong pairing (topologically trivial) 
state gives rise to the dimer phase. In particular the mapping between the Haldane phase and 
Gutwziller wavefunction is exact at the AKLT point K = 1/3. The transition point between the 
two phases determined by the optimized Gutzwiller Projected wavefunction is in good agreement 
with the known result. The effect of Z2 gauge fluctuations above the mean field theory is analyzed. 



I. INTRODUCTION 

Slave boson mean field theory is now accepted as a 
powerful tool in identifying exotic states in strongly cor- 
related electron systemsP*^ At half-filling, the slave bo- 
son approach reduces to a fermionic representation for 
the 5 = 1/2 spins where mean- field theories can be 
built and corresponding trial ground state wavefunctions 
can be constructed through the Gutzwiller projection 
technique.^ The approach has generated large variety of 
wavefunctions used to describe different resonant valence 
bond (RVB) states of frustrated Heisenberg systems in- 
cluding quantum spin liquid statesP^ 

In a recent paperP several of us have generalized the 
fermionic representation to S > 1/2 spin systems and 
have shown that a simple mean-field theory produces re- 
sults which are in agreement with Haldane conjecture 
for the one-dimensional Heisenberg model. A natural 
question is, how about the Gutzwiller projected wave- 
functions obtained from these mean-field states? Are 
they close to the corresponding real ground state wave- 
functions? How about more complicated spin models? 
Here we shall provide a partial answer to these questions 
by studying the Gutzwiller projected wavefunctions ob- 
tained from the mean field st ate s of the 5 = 1 bilinear- 
biquadratic Heisenberg modeP^ 

H=J2[JSi-Sj+K(Si-Sj) 2 }- (!) 

{id) 

where Si are spin operators and j = i + 1 in 1-dimension. 

The 5 = 1 bilinear-biquadratic Heisenberg model has 
attracted much interest. At K = the Haldane conjec- 
ture predicts that the ground state of integer-spin an- 
tiferromagnetic Heisenberg Model (AFHM) is disordered 
with gapped excitations™ Later it was shown by Affleck- 
Kennedy-Lieb-Tasaki (AKLT) that the point K/J =1/3 



is exactly solvable^ and the resulting state is a trans- 
lation invariant valence-bond-solid states.^ The 5=1 
AKLT state together with all states in the so-called Hal- 
dane phase are topologically nontrivial in the sense that 
they cannot be deformed into the trivial S z = trivial 
product state without a phase transition. People have 
been trying to use hidden symmetry breaking,^ or equiv- 
alently, a nonlocal string order to characterize the non- 
trivial order in the Haldane phase. But those character- 
izations are not satisfactory since the Haldane phase is 
separated from the S z = trivial product state even 
when we break all the spin rotation symmetry, in which 
case there is no h idden symmetry breaking and/or non- 
local string order It turns out that the non-trivial 
order in the Haldane phase, called symmetry-protected 
topological order, is described by symmetric local uni- 
tary transformation and the projective representation of 
the symmetry group EES] 

It is now believed that the 5 = 1 antiferromagnetic 
(J > 0) bilinear-biquadratic model in one dimension con- 
tains three phases, the dimer phase at K/J < —1, the 
Haldane phase with — 1 < K/J < 1 and a gapless phase 
at K/J > 1 with antiferro-nematic correlation. 

Dimerized Haldane Critical 
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FIG. 1. The phase diagram of the S = 1 bilinear-biquadratic 
model in the antiferromagnetic section (J > 0). 

The question we shall address in this paper is whether 
the above phase diagram can be reproduced by using 
simple Gutzwiller projected wavefunctions obtained from 
the fermionic mean-field theory. We shall show in the 
following that the optimized projected mean-field wave- 
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functions are very close to the true ground states for 
the ID antiferromagnetic bilinear-biquadratic model in 
the regime K/J < 1. In particular, the optimized pro- 
jected mean field state is the exact ground state at the 
AKLT point K = 1/3. The mean field state is a pair- 
ing state of free fermions, which can be a trivial or 
a non-trivial topological phase, which are classified as 
weak and strong pairing states by their different wind- 
ing numbers.^ The nature of the topological phase of the 
mean field state is found to be important in distinguish- 
ing between Haldane and dimer phases. We find that 
after Gutzwiller projection the weak pairing states be- 
come the Haldane phase whereas the strong pairing states 
become the dimer phase. A long-ranged spin-Peierls or- 
der emerges in the strong pairing states after Gultzwiller 
projection although the spin-Peierls correlation is short 
ranged at mean field level. 

Above results can be understood from the fermionic 
mean field theory when Z^ gauge fluctuations are taken 
into account. We find that the Z% instantons behave 
differently in the weak pairing region and strong pairing 
region. Thus the Haldane phase and the dimer phase can 
also be distinguished by their different effective Z^ gauge 
theories. 

The paper is organized as follows: In section [Hj we re- 
view the fermionic representation for S = 1 spins and in- 
troduce the mean field theory for the bilinear-biquadratic 
model. The Gutzwiller projected mean-field wavefunc- 
tions are introduced in section |III| as trial ground state 
wavefunctions for the bilinear-biquadratic model. Us- 
ing variational Monte-Carlo (VMC) technique^ we find 
that the projected mean-field wavefunctions after opti- 
mization are very close to the true ground states of the 
ID antiferromagnetic bilinear-biquadratic model in the 
regime K/J < 1. We show in particular that the opti- 
mized projected mean field state is the exact ground state 
at the AKLT point K = 1/3. Based on mean-field the- 
ory, we construct in section |IV| the effective low energy 
theories for the bilinear-biquadratic model. The paper 
is concluded in section [V] with some general comments. 
The mapping between the mean-field zero-energy Majo- 
rana end states in the weak-coupling phase and spin-1/2 
end states of open spin chains in the Haldane phase is 
established in the Appendix [A] 



II. FERMIONIC REPRESENTATION AND 
MEAN-FIELD THEORY FOR SPIN S = 1 
BILINEAR-BIQUADRATIC MODEL 

The fermionic representation for S — 1 spins is a gen- 
eralization of the fermionic representation for S = 1/2 
spins. In this representation, three fermionic spinon op- 
erators ci,co,c_i are introduced to represent the S z — 
1, 0, —1 states at each site, and the spin operator is given 
as S a — J2m,n c m^mn c n, where I a is the 3x3 ma- 
trix representation for 5 = 1 angular momentum op- 
erator with a = x,y,z and m,n = 1,0,-1. As in 



usual slave particle method, a particle number constraint 
Ni = E m 4ni c mi = 1 (m = -1, 0, 1) has to be imposed 
on each lattice site i to ensure a one-to-one mapping be- 
tween the spin and fermion states. 

With the single occupation constraint imposed the 
bilinear-biquadratic Heisenberg model (fjj) can be writ- 
ten as an interacting fermion modelpMST 



Xij 



(J -if) A]. Ay]. (2) 



where \ij — J2 m c lni c mj i s the fermion hopping operator 
and Aij — cuC-%j — CojCoj + C-nCxj is the spin singlet 
pairing operator. A mean field theory for this interacting 
fermion model can be obtained by introducing the mean 
field parameters Xij = {x.ij)i &ij — (Ay), and the time 
averaged Lagrangian multiplier A^ to decouple the model 
into a fermion bilinear modeP', 

H mi = -J ^ [XijCmjCmi + h.C.] 

+ (J - *0 £[Ay(ct„4 " 44 + 4i c -ii) + h.c] 



<ij> 



■ Xj&^Cmi + const. 
It is interesting also to introduce the cartesian coor- 

• . 1 / 1 i I 



(3) 



dinate operators, c x = -^=(c_i — c\ 
ci),c z — Co (these operators annihilate the states, \x 

-i\ M\\ in i\ i in\ 



' Z \v) 



x2 



^(|-l) + |l)),|z) = |0),respec- 



tively). In this representation the operators Xij an d Aj 



becomes 



Xij — C xi C xj 



C yi C yj 



A, 



(4) 



and the mean field Hamiltonian reduces to three copies 
of Kitaev's Majorana chain model! 21 ' This representation 
will be used in our later discussion. 

The mean field Hamiltonian H m { can be diagonalized 
by the standard Bogoliubov-de Gennes (B-dG) transfor- 
mation. We shall consider periodic/antiperiodic bound- 
ary condition here. (The case of open-boundary con- 
dition is discussed in the Appendix.) In this case the 
system is translational invariant and Xij = X-, ^ij = 
A, Xi — A become site- independent. The mean-field 
Hamiltonian is diagonal in momentum space, 



H, 



mi 



E 



-[Afc^cL^fe - ^lA-k) + h-c] 



const 



(5) 



rn.k 
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where \k — A — 2Jxcosfc and Afe — — 2i(J — K) A sin k. 
/3 m/ t's are related to c m /fc, c^, fc 's by the Bogoliubov trans- 
formation, 

Pik = u k cik - UfecLi.fe, 

00k = UfeCofc + VfcCo-fc- ( 6 ) 

where = cos %,v k = ?sin%, tan6k = and = 

\/x\ + A^. The mean field dispersion is gapped when 
A^O except at the phase transition point A — 2| J%| = 0. 

The ground state of H m { is the vacuum state of the Bo- 
goliubov particles /3 m fe's. The parameters x an d A are de- 
termined self-consistently in mean-field theory, with A is 
determined by the averaged particle number constraint, 
i.e., 

X=(Xii+x), A=<A« + i>, (JVi)=l. (7) 

where (...) denotes ground state averages. We shall see 
later that the self-consistently determined mean field pa- 
rameters are not optimal in constructing Gutzwiller pro- 
jected wavef unctions. It is more fruitful to treat Xi A 
and A as variational parameters in the trial Hamiltonian 
([3]) that generates a trial mean field ground state Atrial)- 
The Gutzwiller projected mean-field state Pq Atrial) will 
be used as a trial wavefunction for the spin model ([!]). 
The optimal mean-field parameters are determined by 
minimizing the energy of the projected wavefunction 

trial / ■ 

It was shown in Ref . 171 and |2"T1 that the (trial) mean field 
ground state described by Eq.([5]) has nontrivial winding 
number if the mean-field parameters satisfy the condition 

A^O, -2\J X \ <\<2\J X \. (8) 

An important consequence of nontrivial winding num- 
ber is that topologically protected Majorana zero modes 
will exist at the boundaries of an open chain in mean- 
field theory. The mean-field states satisfying Eq. ^ are 
called weak pairing states. The winding number vanishes 
if A 7^ and |A| > 2| Jx\ and these states are called strong 
pairing statesPS We shall show later that the weak pair- 
ing states become the Haldane phase, while the strong 
pairing states become the dimer phase after Gutzwiller 
projection. In later discussion, we will mainly focus on 
the antiferromagnetic interaction case J > 0. To simplify 
notation we shall set J = I in the following. The value 
of J will be defined again only in exceptional cases. 



III. GUTZWILLER PROJECTED 
WAVEFUNCTIONS 

The Gutzwiller Projection for 5 — 1 systems is in prin- 
ciple the same as Gutzwiller Projection for 5=1/2 sys- 
tems. In the mean field ground state wavefunction, the 
particle number constraint is satisfied only on average 



and the purpose of the Gutzwiller projection is to re- 
move all state components with occupancy Ni =/= 1 for 
some sites i, and thus projecting the wavefunction into 
the subspace with exactly one fermion per site. 

There are however a few important technical differ- 
ence between spin-1/2 and spin-1 systems. First of all, 
5 = 1/2 systems are particle- hole symmetric and the 
particle number constraint is invariant under particle- 
hole transformation. As a result, we can always set 
A = in the trial wave functions. This is not the case 
for 5 = 1 models where the particle number constraint 
is not invariant under particle-hole transformation. Con- 
sequently, A 7^ in general and should be treated as a 
parameter determined variationally. Notice that A de- 
termines the topology of the mean field state and the 
corresponding Gutzwiller projected states as we shall see 
in the following. 

The second important difference between spin-1/2 and 
spin-1 systems is that a singlet ground state for 5=1 
systems is composed of configurations with different S z 
distribution (ni,noi n -i) = ( n i^L — 2ni,ni), where n m 
is the number of spins with S z = m in a spin configura- 
tion and L is the number of sites in the spin chain. It 
is clear that ri\ can take any value between and [L/2] 
where [L/2] = L/2 if L is even and [L/2] = (L - l)/2 
if L is odd. On the contrary rii/2 = "--1/2 = L/1 is 
fixed for spin-1/2 systems. Notice that as a result L is 
always even for spin-1/2 singlet states, but can be even 
or odd for spin-1 singlet states. For even L, n is even 
and the Gutzwiller projected wavefunction is a straight- 
forward projection of the mean-field ground state which 
is a paired BCS state. For odd L, no is odd. This means 
there is one 5 Z = fermion mode (cq) remaining un- 
paired in the ground state. The projection is similar to 
the even L case except that we have to keep in mind the 
occupied free fermion mode. The situation for open spin 
chains is further complicated by the existence of Majo- 
rana end modes which is discussed in the Appendix. 

The above condition results in a natural choice of 
boundary condition in constructing the Gutzwiller pro- 
jected wavefunction. To see this we first note that the 
values of allowed fermion momentum fc in periodic and 
anti-periodic boundary conditions are different. They 
are given by k = ( 2M + 1 ) 7r ) under periodic (an- 

tiperiodic) boundary conditions, where M take values 
M = -[^], -[^} + 1, [f ]. The energy spectrum 
Ek is doubly degenerate except at the points k = which 
exists only under periodic boundary condition (for both 
even or odd L), and k = tt which exists under periodic 
boundary condition for even L and under anti-periodic 
boundary condition for odd L. Notice that at these two 
points Afc = — 2iAsinfc vanishes and they are natural 
candidates for constructing a singly occupied fermion 
state. All other momentum fc's are paired at the ground 
state. The energies at these two points are given by 
Eq(tt) = A - (+)2x- Notice that E < and E n > in 
the weak pairing phase (we choose a gauge where x > 0) 
whereas both Eq^ > in the strong pairing phase. 
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TABLE I. Comparison of the energies obtained from VMC (with L — 100) and those from other methods. The points (1, 1)uls, 
(1, g)AKLT, (1, — 1)tb, (0, — l)su(3) ar e exactly solvable. The comparison energy of the rest points are obtained with the 'infinite 
time-evolving block decimation' algorithm (we have set the dimension of the matrices as 48).^ In the last line we list the optimal 
variational parameters (x, A, A) obtained by VMC. 



(J,K) 


(1,1)uls 


(1, ^)aklt 


(I7 0) Heisenberg 


(1,-1]tb 


(1,-2) 


(1,-3) 


(0,-l)suf3) 


(-1,-3) 


(-1,-2) 


comparison 


0.2971^ 


— A — 


-1.4015^ 




-6.7531 


-9.5330 


-2.7969^ 


-7.3518 


-4.5939 


VMC 


0.299#] 
±0.0004 


±7 x 10" 15 


-1.4001 
±0.0004 


-3.9917 
±0.0012 


-6.7372 
±0.0023 


-9.5103 
±0.0034 


-2.7953 
±0.0005 


-7.2901 
±0.0038 


-4.4946 
±0.0028 


(X,A,A) 


(1,0,1) 


(1, |,o) 


(1,0.98,1.78) 


(1,1.11,2.00) 


(1,1.15,2.07) 


(1,1.79,2.22) 


(0,1,0.14) 


(0,1,0.21) 


(0,1,0.12) 



a Due to the 5(7(3) symmetry, the particle number of c x ,c y , c z should be equal. To this end, we have set L = 99. 



We now consider the weak pairing phase. In this case 
a lowest energy state is formed for chains with odd L 
when the k — state is occupied, i.e. periodic boundary 
condition is preferred. On the other hand, for chains with 
even L an anti-periodic boundary condition is preferred 
so that the k = state is not available and a paired BCS 
state is naturally formed. 

Next we consider the strong pairing phase. Suppose 
L is even. Under anti-periodic boundary condition, all 
the fermions are paired in the ground state. Under pe- 
riodic boundary condition, since Eq^ > 0, the unpaired 
fermion modes at k = and k — it are unoccupied. The 
remaining fermions are paired. Both of the two bound- 
ary conditions are allowed. When L is odd, the ground 
state is not a spin singlet and behaves differentlyPS In 
later discussion about the strong pairing phase, we will 
mainly consider the case where L is even. 

The above result suggests that the weak-pairing 
ground state is unique, whereas the strong-pairing 
ground state is doubly degenerate. To see this we note 
that for a closed chain, the existence or absence of tt flux 
through the ring (which corresponds to the anti-periodic 
boundary condition or the periodic boundary condition 
for fermions, respectively) usually result in two degener- 
ate time-reversal invariant fermion states. However, as 
shown above, for a chain with fixed L, we cannot choose 
boundary condition freely for the Gutzwiller projected 
state in the weak-coupling phase, whereas both bound- 
ary conditions are available in the strong-coupling phase. 

In the following we shall report our numerical results 
for the antiferromagnetic bilinear-biquadratic Heisenberg 
model in the parameter range —00 < K < 1, the point 
K = —00 is obtained by taking (J,K) = (0,-1). For 
comparison purpose we have also studied the Gutzwiller 
projected wavefunctions at J < with (J, K) = (—1, —3) 
and (—1,-2). As in S = 1/2 case, the optimization of 
energy of the projected wavefunctions is carried out us- 
ing a variational Monte Carlo method (VMC). We set 
the length of the chain to be L — 100, and has taken 10 6 
MC steps in our numerical work. We note that because 
of the difference between 5 = 1/2 and S — 1 systems as 
discussed above, the VMC for spin-1/2 Gutzwiller pro- 
jection procedure has to be modified for S = 1 models. 
We shall not go into these technical details in this paper. 
We shall first report in section III A our overall numerical 



results and phase diagram which are in good agreements 
with known results as long as the ground state is a spin- 
singlet. Our Gutzwiller projected wavefunction provides 
a good description of the system even around the critical 
point (K = —1) between the Haldane and dimer phases. 
In section [TlI B| wc illustrate analytically that a projected 
BCS state becomes the exact ground state of the model 
at the AKLT point K = 1/3. 



A. Overall results and phase diagram 

Our numerical results of Gutzwiller projection for the 
bilinear-biquadratic Heisenberg model is summarized in 
Table [I] The ground state energy computed from the (op- 
timized) Gutzwiller projected wavefunction is compared 
with the exact result or results obtained from the 'infi- 
nite time-evolving block decimation' algorithm.^ We see 
that agreement in energy is better than 0.8% in the whole 
range —00 < K < 1(J = 1), where the ground state is a 
spin singlet. The optimized parameters given in the ta- 
ble suggests that the system is in the weak pairing phase 
when — 1 < K < 1 and is the strong pairing phase other- 
wise. The critical point at K ~ — 1 will be studied more 
carefully in the following. To understand the nature of 
the weak- and strong- pairing phases we first study the 
spin-spin or dimer-dimer correlation functions at K = 
and K = -2. 

At K = the self-consistent mean field solution gives 
X = 0.5671, A/x = 1.3447, A/x = 1.3594. The energy of 
corresponding Gutzwiller projected wavefunction is E g — 
— 1.3984 ± 0.0004 per site, which is already quite close 
to the known ground state energy E g = — 1.40ld^2£l 
with a difference ~ 0.2%. The correlation length deter- 
mined from a fitting to the spin-spin correlation function 
is roughly 3.25 units of lattice constant, which is smaller 
than the value 6.03 given in literaturepSH23 

Our result can be further improved by optimizing the 
parameters x, A, A. The optimal parameters we ob- 
tain are x = 1,A = 0.9767, A = 1.7810 (here we nor- 
malize x = !)• The energy of the projected state is 
E g = -1.4001 ± 0.0004 which is further improved by 
0.1%. The spin-spin correlation function is plotted in 
Fig. [2J The correlation (SfS* +r ) matches very well with 
(SfSf, r ), indicating the rotational invariance of the pro- 
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FIG. 2. (Color online) The spin-spin correlation function 
C a (r) — (SfSj+r) in the optimized projected wavefunction 
1 1,0.9767, 1.7810) at K = 0. The staggered sign (-l) r shows 
the short-range antiferromagnetic order of the ground state. 
The error bar of C x (r) is smaller then that of C 2 (r), because 
in calculating the former, more spin configurations are in- 
volved. The dashed lines are exponential fitting of the data. 
The fitting shows that the correlation length is about 5.92 
units of lattice constant. 



jected wavefunction. The correlation length determined 
from the optimized wavefunction is 5.92 lattice constants, 
which is very close to the accepted value 6.02. 

A trademark for the Haldane phase is the existence 
of spin-1/2 end states. Indeed, end Majorana fermion 
states are observed to exist in the weak-pairing phase of 
the fermionic mean-field theory.^ The question is whether 
these Majorana end states become spin-1/2 end states 
after Gutzwiller projection. This question is discussed in 
the Appendix where we show how the Majorana fermion 
end states turn into spin-1/2 end states after Gutzwiller 
projection. 

Next we consider K = —1. In this case, the mean 
field solution has x = 0, A = 0.9203, A = 0.9110, and 
the energy of the projected mean-field state is E g — 
—6.6691 ±0.0023. The energy can be lowered by optimiz- 
ing the mean field parameters. The optimal parameters 
found from the VMC are \ = 1, A = 1.1532, A = 2.0661 
with energy E g = -6. 7372± 0.0023. The Gutzwiller pro- 
jected wavefunction is obviously translationally invariant 
and does not explicitly break the translation symmetry. 
To see that the state described a dimer state, we com- 
pute the spin-Peierls correlation function. The result is 
shown in Fig. [3] We note that the spin-Peierls correlation 
is clearly short-ranged in the Heisenberg model(A~ = 0) 
but is long-ranged in the K = —2 model. Notice that 
the "weak-pairing/Haldane" and "strong pairing/ dimer" 
mapping is in agreement with the ground state degener- 
acy we deduced in last section. It is remarkable that 
the spin-Peierls correlation becomes long-ranged only af- 
ter the Gultzwiller projection and is short-ranged before 
projection. 

Now we examine the phase transition point between 



FIG. 3. (Color online) Spin-Peierls correlation functions at 
K — —2. The blue circled line shows the Spin-Peierls correla- 
tion function (P? P* +r ) of the projected self-consistent mean 
field ground state P G \ip(x = 0, A = 0.9203, A = 0.9110)} mf 
with K — —2. Here the spin-Peierls order is defined as 
Pi — SiSi + i — Sf +1 S , f +2 . The red squared line shows the 
Spin-Peierls correlation of the projected optimal trial mean 
field state ffc#(x = 1, A = 1.1532, A = 2.0661)) tr iai with 
K = —2. For comparison, the black crossed line is the spin- 
Peiels correlation of the projected optimal trial mean field 
state of the Heisenberg model (K = 0). 



the Haldane and dimer phases. Fig. [4] shows the spin 
Peierls correlation at distance L/2 as a function of K 
for the projected optimal trial wavefunctions near the 
phase transition point. Within numerical error our re- 
sults show that A— 2\ = and the spin-Peierls correlation 
vanishes at the point K « — 1. The spontaneous break- 
ing of translation symmetry indicates that the transition 
is of second order, consistent with exact solution at the 
Takhatajan-Babujian (TB) point K = -P We would 
like to point out that the transition point determined by 
the self-consistent mean field theory is K = —0.33, which 
is far away from the exact result. 



B. AKLT point: the projected wavefunction as 
exact ground state 

The validity of the Gutzwiller projected wavefunction 
approach for 5 = 1 spinchains is further supported by an 
exact result at the AKLT pointP' (K = 1/3), where we 
find that the AKLT state can be exactly represented as a 
Gutzwiller projected BCS wavefunction with parameters 
X = 1,A = 3/2,A = 0. 

It is convenient to adopt the cartesian coordinate in 
the discussion. Firstly we consider the closed bound- 
ary condition. Since the MF Hamiltonian H m f in (|3|,([4]) 
contains three identical copies for the flavors c Xl c v ,c Zl 
we may concentrate on a single flavor H a and define 



l 



ilai + i-Tcxi), where 7^ and 7^ are Majo- 
rana operators satisfying {7^,7^} = S^Sapdij where 
a,b = l,r and a,/3 = x, y, z. Setting x — (1 — K)A = 1, 
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FIG. 4. (Color online) Projected optimal trial mean field 
states near the phase transition point. Here we set J = 1 
and x = 1- The spin-Peierls correlation qualitatively shows 
the expected result that A — 2\ have differen sign in different 
phases. When A — 2\ > 0, the projected wavefunctions have 
finite spin-Peierls correlation, indicating they are dimerized. 
Otherwise, when A — 2\ < 0, the spin-Peierls correlation ap- 
proaches zero, indicating the states are not dimerized. The 
point where A — 2\ = almost overlaps with the point where 
(P* P z ,l) — ¥ 0. This result verifies our conclusion that the 
topology of the states distinguishes different phases. 



A = and adopting close boundary condition \li = 
(1 — K)Alx = (— l) ri (rj = for periodic boundary condi- 
tion and rj = 1 for antiperiodic boundary condition), the 
mean-field Hamiltonian H a can be mapped into Kitaev's 
Majorana chainj^ 

L-l 

H a = (-l)"(-2z 7 ^ 7 ii) + Yi-ZKilL+i), (9) 



where we have dropped an unimportant constant. Notice 
that every term in ^ has eigenvalues ± 1 and is commut- 
ing with all other terms. Defining the fermion parity^ 

Pa = e -£^ = Y[{1 2clc ai ) = 11(2*7^), 

i i 

it is easy to see that the fermion parity of the ground 
state of the mean-field Hamiltonian is given by 

p a = -2 l7 r aL j l al n( 2 *£iVai +l ) - (io) 



and is even(odd) under anti-periodic (periodic) boundary 
condition. The total fermi parity is the product of the 
three flavors P fermi = P x P y P z = {P x ) 3 = P x . Notice 
that the periodic boundary condition is ruled out by the 
particle number constraint for even L. This is a general 
property of the BCS Hamiltonian in weak-coupling phase 
as we have discussed before. 

Since all terms in ^ are commuting, the ground state 
of one term in the Hamiltonian ^ provides a supporting 
Hilbert space for the reduced density matrix constructed 



for the whole ground state. The ground state of H a for 
two neighboring sites i and i + 1 is two-fold degenerate, 

\<Mi = ( c L + c L+i)l vac >> 
\<f> a }2 = ( 1 + c ii c L+x)W ac )- 

Since there are three flavors a = x,y, z, the ground state 
of the two sites is a product state \<j>) — l^)!^)!^), 
and are 8-fold degenerate. It is easy to see by direct 
computation that the Gutzwiller projection kills half of 
these states, and the surviving four states are 



(J J 

(J J 



xi zi-\-l 



^y&xi+l 

L 'zi L 'xi+ 
t J 



+ r 1 c 



t J 



zi zi+1 



vac) 



vac) 



i)|vac), 



cue. 



yi+li 



The first one is a spin singlet, and the remaining three 
form a (S = 1) triplet. The absence of spin-2 states for 
every two neighboring sites is a fingerprint of the spin-1 
AKLT stateFCQ Thus, we prove that the Gutzwiller pro- 
jected trial state with x — 1, A = 3/2, A = is equivalent 
to the spin-1 AKLT state. 

We note that the above proof can be extended straight- 
forwardly to the SO(n) symmetric AKLT models.^ 
When n is odd, the n Majorana fermions at the edge form 
the irreducible spinor representation of SO(n) group. So, 
after Gutzwiller projection, the ground state of n-copies 
of Kitaev's Majorana chain model exactly describes the 
ground state of the S'0(rt)-AKLT model. When n is even, 
the ground state is dimerized since the spinor representa- 
tion of SO(n) is reducible. The equivalence between the 
ground state of the S'0(n)-AKLT model and the pro- 
jected ground state s of n-copies of Kitaev's Majorana 
chain remains validPEH In that case, the n Majorana 
fermions at the edge form a direct sum of two versions of 
irreducible SO(n) spinor representations. 



IV. EFFECTIVE LOW ENERGY THEORY 

The success of the Gutzwiller projected wavefunction 
in describing the ground state properties of the Haldane 
and dimer phases suggests that the low energy properties 
of these phases may be well described by effective field 
theories constructed from the corresponding mean-field 
states. We provide two approaches in this section. The 
first one is a Z 2 gauge field description by integrating out 
the fermions, and the other one is an effective (Majorana) 
fermionic field theory. 

The mean-field theory provides a correct description 
of the low energy properties of the Haldane and dimer 
phases when Z 2 gauge fluctuations or Z 2 instantons are 
taken into account. Since the mean field state is a fermion 
paired state with finite gap, the resulting low energy ef- 
fective theory is Z 2 gauge theory. Usually, the Z 2 gauge 
theory in (1+1)D is always confined since the Z 2 instan- 
tons (ie the Z 2 vortices in (1+1)D discrete space-time) 
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have a finite action. However, the low energy effective 
Z2 gauge theory obtained from our models has different 
dynamical properties. 

In the weak pairing phase, a Z2 instanton gives rise 
to a fermionic zero modeP^ As a result, the action of 
separating two instantons is proportional to —t/£, where 
t is the time-distance between the instantons and l/£ 
is the excitation gapP^ This means that the Z2 instan- 
tons are confined and consequently the Z 2 gauge theory 
is deconfined (this situation still holds if the Hamiltonian 
contains a dimerized interaction) . The Z2 vortex changes 
the fermion parity of the mean field ground state. Ow- 
ing to the particle number constraint, a permitted in- 
stanton operator should be a composition of a Z2 vortex 
and a spinon operator. Thus, an instanton carries 7r- 
momentum and spin-1. An example of such instanton 
operator is given as 

<p = ]T(-iys 4 = s w . 

i 

Above instanton operator creates a magnon with momen- 
tum 7T. 

In the strong pairing phase, due to the absence of 
fermion zero modes, the Z2 vortices in (1+1)D space-time 
have a finite action and consequently have a finite density. 
However, since the Z2 instanton carry tt crystal momen- 
tum, there will be an extra phase factor (—1) associated 
with it. When sum over the contribution of instantons at 
all spacial positions, the phase factors will cause cancela- 
tion. Consequently, the effect of instantons is suppressed, 
and the Z2 gauge theory is still deconfinedPS On the 
other hand, if the Hamiltonian have a translation sym- 
metry breaking term (such as a dimerized interaction), 
then the action of the instantons will not be canceled and 
will confine the Z2 gauge filed. (This is a remarkable dif- 
ference between the strong pairing phase and the weak 
pairing phase.) Since a Z2 instanton carries 7r momen- 
tum and zero spin, we can give an example of such an 
operator 



n+i 



(-1) 



i+i 



The finite action of the instantons indicates that the 
ground state is degenerate and has a finite spin-Peierls 
order. 

At the transition point, the mean field dispersion is 
gapless at k = due to A— 2\ = 0, and the low energy ex- 
citations are consist of three species of Majorana fermions 
with energy E k w 2(1 - K)A\k\. Notice that three Ma- 
jorana fermions form a spin-1/2 objectP^ Consequently, 
after Gutzwiller projection, the elementary excitations 
carry spin-1/2. Notice that the number of excited Ma- 
jorana fermions must be even in a physical state, so the 
spin-1/2 excitation must appear in pairs. This physi- 
cal picture arising from mean-field theory agrees well (in 
long wavelength limit) with result coming from the Bethe 
ansatz solution of the TB modelPS 



Since the effective Z2 gauge fields are deconfined in 
both the Haldane and the dimer phases of model ([l}, it 
will be a good approximation to ignore the gauge field 
and consider the fermion theory only. Tsvelik proposed 
an effective Majorana field theory to describe the low- 
energy physics close to the TB poinlP^ 



K 



eft 



2~^ 



- l l a d xl l a ) - im^lall (11) 



where 7^ and 7^ are right and left moving Majorana 
fermions. Marginal terms (four-fermion interactions) are 
neglected. This theory describes the Haldane phase for 
m < and the dimer phase for m > 0. Thus, the quan- 
tum criticality at the TB point belongs to Ising univer- 
sality class. Our mean field Hamiltonian ( J5| i n Majo- 
rana fermion representation is the same as (111 in long 



wave length limit (strictly speaking, one can only com- 
pare the mean field theory with the effective field the- 
ory after renormalization) . The fermion mass and the 
velocity are related to the mean field parameters up to 
renormalization factors 

mocA-2J|x|, v oc 2( J — K)A. 



Strictly speaking, the effective theory |TT| ) is only valid 
when the spin Hamiltonian is translationally invariant. 
Otherwise, the Z2 gauge field is confined in the dimer 



phase, and (111 will not describe the low energy behaviors 



near the transition point correctly. 



V. DISCUSSION AND CONCLUSION 

We give a few comments about the regimes K > 1 (J — 
1) and J < in the following. 

In the region K > 1 the pairing term in ^ becomes 
irrelevant and consequently A = in our mean-field the- 
ory. In this case, the (trial) mean field ground state has 
a 1/3-fillcd fcrmi sea, whose fermi points are located at 
k-p = ±7r/3. Physically, for K > 1, the marginally irrele- 
vant instability of the model ([I]) is the antiferro-nematic 
order PES! Notice that there should be no true long-ranged 
antiferro-nematic order in the ground state in one dimen- 
sion. Nevertheless antiferro-nematic correlation is not 
included in the mean field ansatz we propose here and 
we do not expect the corresponding Gutzwiller projected 
wavefunctions will describe the ground state well in this 
regime. One need to introduce new mean field parame- 
ters which is beyond the scope of the present paper. 

The situation may be different at the special point K = 
1 which corresponds to the integrable SU(3) ULS model. 
Bethe ansatz solution at this special point indicates that 
excitations above the SU(3) singlet ground state are 
gapless^ and are described by a SU(3)i Wess-Zumino- 
Novikov-Witten (WZNW) model with marginally irrel- 
evant perturbations!^ The physics of the SU(3) ULS 
model can be obtained from the Gutzwiller projected 
wavefunction. The pairing term in ([3| vanishes at J = 
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K = 1 and the MF Hamiltonian becomes a free fermion 
model. Moreover A = \ since the fermion bands are 
1/3 filled with fermi points at k = ±7r/3. In this case 
the variational wavefunction |Whs) — -PgIV'(X) Oj x))mf is 
an SU(3) singlet and is the exact ground state of the 
SU(3) Haldane-Shastry model with inverse-square inter- 
actions.^ Since the effective field theory of the SU(3) 
Haldane-Shastry chain is just an unperturbed version of 
the SU(3)i WZNW modelp^ which shares similar low- 
energy physics with the ULS model, we expect that the 
Gutzwillcr projected fermi sea state is a good variational 
wavefunction for the K — 1 antiferromagnctic bilinear 
biquadratic model (see also Tab. [I]). 

Lastly we consider the J < regime. There are two 
phases. The ferromagnetic phase is located at J < K < 
oo. This phase is again beyond our present mean field 
theory which is designed to describe spin-singlet states. 
The remaining part K < J also belongs to the dimer 
phase and is described by the projected BCS state in 
strong-pairing regime. The difference between this dimer 
phase and the dimer phase at J > is that the hopping 
term in ([2]) becomes irrelevant and x vanishes when J < 
(see TabTllfor three examples). 

To conclude, we introduce in this paper Gutzwiller pro- 
jection for the mean field states obtained from a fermionic 
mean-field theory for S = 1 systems. The method is ap- 
plied to study the one-dimensional bilinear-biquadratic 
Heisenberg model. We find that the topology of the mean 
field state determines the character of the Gutzwiller pro- 
jected wavefunction. The projected weak pairing states 
belong to the Haldane phase while the projected strong 
pairing states belong to the dimer phase. This result 
is consistent with the Z-i gauge theory above the (trial) 
mean field ground state. Our theory agrees well with the 
Majorana effective field theory at the TB critical point 
and the method can be generalized to higher dimensions 
or other spin models. 
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Appendix A: End states in open Heisenberg chains 

It was shown in Ref.[7]that the mean field ground state 
in the fermionic mean-field theory for S = I Heisenberg 
model is topologically nontrivial. As a result, each of the 



three flavors of fermions c x ,c y , c z have one Majorana zero 
mode located at each end of an open chain with wave- 
function exponentially decay away from the chain end. 
(We adopt the cartesian coordinate here). It was shown 
in Ref. [33] that three Majorana fermions together may 
represent a spin-1/2 object via S a — je 0135 ^ 13 ^ 6 , where 
a, /3, 5 = x, y, z and 7" is a Majorana fermion operator 
corresponding to the flavor c a (see also section IIIB). 
This suggests that the mean field Majorana zero modes 
may form spi n-1/2 sta tes at each end, in agreement with 
known resultsP- 5 * 26 * 40 ! In this appendix we show how the 
Majorana zero modes become the spin-1/2 edge states af- 
ter Gutzwiller projection. First we demonstrate numeri- 
cally the existence of Majorana end modes in mean-field 
theory for open spin chain. 



LU 




200 



FIG. 5. (Color online) The mean field excitation spectrum Ei 
of each flavor for an open chain. Due to particle-hole sym- 
metry, the 200 modes correspond to L — 100 fermion modes. 
The two Majorana zero modes manifest themselves. The en- 
ergy of these two Majorana modes exponentially approaches 
to zero with increasing L. 

The existence of majorana fermion edge mode in mean 
field theory can be seen from the mean-field dispersion Ek 
for open chains (Fig. [5]), where the existence of two zero- 
energy modes (for each favor) are clear. These two Majo- 
rana modes, noted as 7^ and 7^ can be combined into a 
zero-energy complex fermion mode c° nd = 7^ + i"f l a which 
can either be occupied or unoccupied. Correspondingly 
the fermion number of the flavor c a fermion can be either 
even or odd (fermion parity) in the ground stated We 
note that strictly speaking, under open boundary condi- 
tion, the mean field parameters Xij j ^ij an d Aj vary from 
site to site and should be determined self-consistentlyP^ 
We assume for simplicity that they take uniform values 
determined by closed spin chain. 

For the three flavors of fermions, there are a total of 
six Majorana zero modes, resulting in 8-fold degenerate 
mean field ground states. Half of these states have odd 
total number of fermions, and half have even. For chains 
with fixed length L, only half of them survives for the 
same reason as discussed in section III. 

The remaining four states have different fermion par- 
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ity distributions (P x , P y ,P z ). One of them is (even, even, 
even). In this case, the three flavors have equal weight in 
the projected state, corresponding to a spin-singlct. The 
other three are (even, odd, odd), (odd, even, odd) or 
(odd, odd, even). These three states form a spin-triplet 
because they can be transformed from one to another by 
a global spin rotation. These four states remains degen- 
erate (in the L — ¥ oo limit) after Gutzwiller projection, 
suggesting that the majorana fermion states become sta- 
ble S = 1/2 end spin states after projection. 



FIG. 7. (Color online) Projected spin-triplet ground 
state with fermion parity distribution (N x , N y , N z ) — (odd, 
odd, even). The spin-spin correlation function of 
{Sf^ N Si^ N+1 ^ L } w -1/4 when TV > 15 shows that the 
edge states carry spin-1/2. 




FIG. 6. (Color online) Projected singlet ground state. The 
spin-spin correlation function of {Sl^ N Sf J _ N+1 ^ L ) ~ —1/4 
when N > 15 shows that the edge states carry spin-1/2, here 

Sl^N = Ei=l &i an d S'£_ J v+l^L ~ Ej=l &L+l-j- 




To confirm the above picture we calculate the cor- 
relation function (S±^ N Sf j _ N+1 ^ L ) , where Sf^jy = 

Eili S f and S l-n+i^l = T,f=i S l+i-j measures the 
total spin S z accommodated from sites i = 1 to JV and 
from i = L — N + 1 to L, respectively. Note that 



(Sf. 



(Sf^. N Sf j _ N+1 ^ L ) in the singlet state (even, even, even). 
The result is shown in Fig. [6] where we find that the cor- 
relation function approaches —\ when N > 15, indicat- 
ing that the magnitude of effective end spins (S^ N and 



and (S Z L _ 



0. We first compute 



5' 



L+l-JV- 



is exactly 1/2. 



The same analysis can also be applied to the 
triplet states. We consider in Fig. [7] the correlation 
(^i^n^l-n+i^l) m tne (odd, odd, even) state. It is 
clear that (S z ^ N Sl_ N+1 ^ L ) m -1/4 when N > 15, 
confirming the existence of spin-1/2 effective end spin in 
the weak-pairing phase. 
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